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Catalysis of nuclear reactions by metastable charged particles X~ opens the possibility for pri- 
mordial production of elements with A > 7. We calculate the abundance of ^Be, where synthesis 
is mediated by the formation of (*BeX~) bound states, finding a dramatic enhancement over the 
standard BEN prediction: ^Be/^H ~ 10~^^ x (Yx/10~^). Thus observations of ^Be abundances at 
low metallicity offers a uniquely sensitive probe of many particle physics models that predict X~ , 
including variants of supersymmetric models. Comparing the catalytically-enhanced abundances of 
primordial ^Li and ®Be, we find the relation ^Be/^Li — (2 — 5) x 10~^ that holds over a wide range 
of X~ abundances and lifetimes. 



The first model of Big Bang Nucleosynthesis [H (BBN) 
made an ambitious attempt at explaining all elemental 
abundances as a result of successive primordial neutron 
capture. Two of the most important reasons why that 
theory did not work, and only very light nuclei can be 
generated in the Big Bang, are the nuclear A ~ 5 and 
A — 8 divides, or the absence of stable nuclear isotopes 
with these mass numbers. After five decades of progress 
in nuclear physics, astrophysics and cosmology, we now 
have a very successful framework of Standard Big Bang 
Nucleosynthesis (SBBN) that makes predictions for el- 
emental abundances of light elements, H, D, He and 
Li, as functions of only one free parameter, the ratio of 
baryon to photon number densities 0. With an addi- 
tional CMB-derived [3] input value of 77^ = 6 x lO^^", 
the comparison of observed amounts of deuterium, he- 
lium and lithium serve as a stringent test of Big Bang 
cosmology and particle physics, and the importance of 
these tests is paramount for many theories seeking to ex- 
tend the Standard Model 0. 

Last year it was realized that BBN is sensitive to a 
wider class of particle physics models than was previously 
thought through the phenomenon of catalyzed primordial 
nuclear fusion or CBBN. In particular it was shown 
that the presence of metastable heavy negatively charged 
particles, that unavoidably will form bound states with 
nuclei H, 0], leads to the catalysis of ^Li production 
via (^HeX") -f D *^Li + X" 3, with a rate exceeding 
the SBBN rate by many orders of magnitude. Moreover, 
the earlier presence of these particles at T ~ 30 keV 
leads to a moderate reduction of ''Li abundance [3 , thus 
tantalizingly reproducing the observational pattern of ^Li 
and ^Li [8|. Since then many groups have incorporated 
the catalyzed reactions into their calculations , and the 
first dedicated three-body potential model calculation of 
the catalyzed rate was performed in [loj . 

The catalysis of ^Li can be viewed as an effective 
"bridging" of the A = 5 divide, and in this Letter we 
show that = 8 is also bridged resulting in an enor- 
mous enhancement of the primordial ^Be abundances 
over the SBBN value, ^Be/^H < lO-^^ thus ef- 



fectively incorporating ^Be into the BBN "family" of 
light nuclei. As was pointed out in the original papers 
U 01, the path to A > 8 nuclei is controlled by the 
{^BcX~) bound state, to get to which X~ has to go 
through a "double bottleneck" of successive a-captures: 
X-^{'^RcX-)^{^BcX-). Once the (^BeX") bound 
state is formed, it participates in neutron capture, 

CBBN: {^BeX-)+n~> ^Be + X"; = 0.26MeV, (1) 

that is catalyzed relative to neutron capture on ^Be, 

SBBN : *Be + n ^ ^Be -I- 7; Qs = 1.665MeV. (2) 

The latter is, of course, never important for BBN because 
free ^Be lives for under a femtosecond. According to the 
general scaling of quantum mechanics [3], the cross sec- 
tion for ([1]) is enhanced relative to (l2|) by a large factor 
~ c^/{uja)^, where oj = Qs/fi is the photon frequency in 
([2]), a is the characteristic size of the {^BeX~) system, 
and c is the speed of light {c = h — 1 hereafter, and 
a — e'^). In the rest of this Letter we analyze the rate 
for reaction ([T]) in detail, calculate the rates for (*BeX~) 
formation, and incorporating them in the reaction net- 
work, calculate the resulting ^Li, ^Be synthesis at T ~ 8 
keV as a function of initial A" "abundance per baryon Yx, 
and lifetime tx- In this work we use the following input 
values for the r.m.s. charge radii of "'He, '^Be, and ^Be 
and calculate the binding energies of (NA") in the limit 
of very heavy A "using a Gaussian charge distribution: 

rSHc = l-67; rfBo = 2.50; r^^Bc = 2.50 fm (3) 

^'ho = 347; El^^ - 1408; EI^^ - 1477 keV. (4) 

While ("^HeA") binding energy has an uncertainty under 
a few keV, the (correlated) uncertainty in Eg^^ and Eg^^ 
can be as large as 50 keV, if e.g. the charge distribution 
is varied from Gaussian to square. More precise values of 
binding energies are possible through theoretical and ex- 
perimental advances in determining the electromagnetic 
form factors of Hght nuclei JJi] . 

Catalyzed neutron capture. It is well known that the 
^Be -|- n system has a low-energy neutron s-wave reso- 
nance {En — 70keV that corresponds to a 1/2+ 
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FIG. 1: Isobar diagram for the *Be + n + X~ system. The 
neutron capture occurs via a (^Be^^X") resonance with subse- 
quent decay to a p-wave continuum of X~ and ^Be in its ground 
state (|") releasing Q = 257 keV. The arrows show the path for 
the CBBN reaction that dominates over 7-emission to the ground 
state of (^BeX~) by four orders of magnitude. 



excited state of ®Be. Through the photodisintegration of 
^Be, the energy and the width of this excited state are 
determined very accurately [l3| . 



dimension [distance] Rio is the "atomic" Is radial 
wave function of the decaying (^BeX~) system, Rpi is 
the final state I — 1 wave function with p = y^2m9^^Q, 
normalized to the p/27r scale , that in large r limit be- 
comes the wave function for the Coulomb problem, /(r) 
represents a "form factor" that we have to assign to the 
interaction of a nuclear dipole operator with the electric 
field created by X~, such that f{r) 1 if the orbit 
of ^Be were far from X^. For our calculation, we take 
/(r) = Att /q x'^dxpc{x)/e, so that for a constant charge 
density it scales as / ~ inside the nuclear radius and 
f — I outside. Explicit numerical calculation gives 
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(8) 



We note the resulting length scale in ([8]) is six times larger 
than the effective size of the nuclear dipole in ([6]) giving 
some a posteriori justification to our procedure. Combin- 
ing (l6|), (O, and ([8]) we arrive at the following estimate 
for the decay width, 



Ei+ = 1735±3keV; r„(£;„) ~ 2 
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keV ' 



r^, = 0.57eV. 



(5) 
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Notice the \/E behavior of the width for the near- which constitutes four orders of magnitude enhancement 



threshold resonance 1J| . In Fig. [T] we show the isobar 
diagram for the *Be + n + X~ system, that holds the 
key for determining the rate of reaction ([1]). As is ev- 
ident from this diagram, the neutron resonance occurs 
right at the threshold, E^^"" ~ ± 30keV, with the esti- 
mated error related to differences in charge distribution 
for ^Be and *Be. Fortunately, the neutron width in ([5]) 
is larger than possible shifts of resonance energy, which 
makes this uncertainty tolerable for the calculation of 
the rate (even if the actual resonance is slightly below 
the threshold). Furthermore, there is no reason to as- 
sign a different width for the neutron resonance in the 
(^BeX~) -I- n system relative to ®Be -I- n and we adopt 
^n{En) from Eq. ([5]). An accurate calculation of the 
decay width for the (^Be^ X^) Be| + X^ process 
requires a dedicated nuclear many-body calculation that 
we will not pursue in this paper. Instead, we choose to 
re-process the existing experimental information on F^ 
by extracting the strength of the reduced matrix element 
for the nuclear El transition. 



over T^. This estimate looks natural, perhaps on the 
lower side, given that Tout contains no small parameters 
and the energy level spacing in this system is on the order 
of 1.5 MeV. Defining rtot(^^n) = Tout + r„(£;„), we cal- 
culate the cross section of catalyzed photonless neutron 
capture ([l]) using the Breit-Wigner formula. 
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g is the spin factor, and g — 1 for llj. The result is a 
factor of '--^ 5 lower than the unitarity limit at £' '--^ lOkeV. 
Thus at temperatures T/IO^K = Tg = 0.1 the reaction 
rate is calculated to be 
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and connect it with an approximate expression for the 
(^Be^~''X~) p-wave decay width, 

Fout ^ -x|(l/2+||d/e||3/2-)px42; (7) 
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r^drxRw{r)^Rpi{r). 



In Eqs. (O and ([7]) ?; is the velocity of the outgoing ^Be 
in the catalyzed reaction, Ir is the radial integral with 



which is large but not larger than e.g. the rate of neutron 
capture on ''Be. Since the actual position of E^^^ has an 
uncertainty of ±30keV we do not have sufficient precision 
to obtain the variation of the rate with temperature. We 
did check, however, that variation in the resonant energy 
by 30 keV introduces only moderate shifts of a fac- 
tor of ~ 2. Variations in Fout directly affect the rate. 
A moderate increase of Fout by a factor of 4 would in- 
crease the cross section to near the unitarity bound while 
a decrease of Fout will obviously have an opposite effect. 
Fortunately, dedicated ab-initio calculations of F„(ii^„) 
and Fout are feasible in state-of-the-art nuclear physics 
[l5| . with the clear potential of improving the accuracy 
of estimate (fTTj). 

Formation of (®BeX~). Before neutrons can undergo 
capture and form ^Be, first the bound states (^IIeX~) 
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TABLE L Properties of (*BeX~) bound states with respect 
to nl quantum numbers. Binding energies are shown relative 
to the *Be+X~ continuum. Resonant energies are given for 
the (*HeX-)+''He system. 

and then {^BcX^} have to form, so that the path to 
^Be is guarded by the double bottleneck. There are 
two main paths to {^BeX~): linear in Yx due to the 
a-reaction on ("^HeX"); and quadratic in Yx due to for- 
mation of the neutral molecules {^lleX~ X^) that react 
with helium via a strong Coulomb-unsuppressed process 
{^neX-X-) +4 He ^ (^Be^-) + X- . Both channels 
first require the formation of (^HeX^) that cannot form 
in any significant amounts above 10 keV due to photodis- 
sociation. We consider these mechanisms in turn. 

In the interesting range of temperatures, the radiative 
fusion of ''He on (^HeX~) proceeds via the formation of 
the resonant atomic states with n = 3, similar to what 
is found in [t^. ®Be has an energy excess of 92keV rela- 
tive to two alpha particles. This leaves only the n = 1,2 
levels below the threshold of ('*HeX^)-f'*He continuum. 
The energy levels of the relevant {^BeX~) excitations, 
resonant energies for incoming a-particles and the elec- 
tromagnetic decay widths are given in Table 1. The key 
observation that facilitates the whole treatment is the 
validity of the narrow resonance approximation. By an 
appropriate rescaling of the Gamow factor for free ®Be 
decay, it is easy to see that at relevant energies the en- 
trance widths Fill to 31 states are ^ O(keV) and thus sat- 
isfy an important condition, T-y <^ Fin <C T, which makes 
the cross section independent of Fin and totally deter- 
mined by Ercs and F-^. Choosing 5 = 2/ -I- 1 in the Breit- 
Wigner formula, and retaining the contributions from the 
two most important resonances, 3d and 3p, we derive the 
thermal rate (av) for the C^HeX") -h^'He (^BeX") -^7 
reaction {Q ~ 969keV) to be 

^q5j.-3/2 (0.95 exp[-1.02/T9] +0.66 exp[-1.32/r9]) . 

(12) 

The non-resonant part of the rate can be somewhat en- 
hanced compared to a typical (a, 7) reaction due to the 
subthreshold 21 resonances, but at T ~ lOkeV it is totally 
negligible in comparison with (fTS]) . At these tempera- 
tures and with this rate the capture of ''He on C^BcX^) 
is about two orders of magnitude slower than the Hubble 
rate and rapidly dropping with T. This is also important, 
as it shows that only a relatively small but non-negligible 



fraction of (^HeX~) will be converted to {^BeX~). 

The rate depends very sensitively on the 3p and 3d 
energy levels, but subtleties of the charge distribution 
in *Be make little difference for their energies. We find 
it remarkable that one can calculate the abundance of 
{^BeX^) virtually free of major nuclear physics uncer- 
tainties. In fact, the main correction to p2)) comes from 
the msBe/wx-suppressed contribution to resonant ener- 
gies that we ignore in this paper, but can account for 
very easily. For X^as light as 100 GeV, these correc- 
tions amount to a 12 keV upward shift of Ey-es, resulting 
in a factor of a few suppression in {^BeX~) abundances, 
but quickly become negligible for heavier mx- 

The molecular mechanism of forming (^BeX~) is com- 
pletely different, as it is regulated by Coulomb unsup- 
pressed processes. Due to the Y^ scaling, this mecha- 
nism is of secondary importance because Yx is rather 
tightly constrained by ^Li overproduction We calcu- 
late the rate of molecular formation (^HeX~) + X~ ^ 
i^BeX^X-) +^ Q ~ 320keV, in the spirit of Kramers 
and ter Haar [l6|, treating the "nuclear" motion of 
semiclassically, and "electron" motion of ^He quantum- 
mechanically, which is known to give a good approxima- 
tion to a full quantum mechanical treatment. At tem- 
peratures Tg <ti Q one finds 

M„,oi = Sn'/'T-'/^J r^dr x T^{r)\V{r)\^'\ (13) 

where r is the distance between two X^ particles, V{r) 
is the potential energy of {^BeX^)-X^ interaction, and 
F^(r) is the probability per time for a quantum jump of 
"^He from the atomic to the molecular state with emission 
of a photon. Using the variationally determined molecu- 
lar wave functions, we calculate F^(r), V{r) and find the 
following estimate for the rate of molecular formation 
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Note that the molecular rate is significantly smaller than 
the "atomic" {^BeX~) recombination rate, 8x lO^Tg"^^^, 
but is not suppressed by rrniie/mx, which is a direct con- 
sequence of the Coulomb attraction in the initial state. 

Synthesis of ^Be. Rates (HH), and ^ enable 

us to calculate the ^Be freezeout abundance numerically. 
Before we do that, we would like to mention that in the 
hypothetical limit of both Yx and tx being large the 
production of ^Be will be neutron-supply-limited, and all 
neutrons produced in the DD and DT fusion below Tg = 
0.1 may end up captured by {^BeX~) before they decay, 
leading to ^Be[large Yx,tx]'^ 0(10~®). Given that at 
lowest metallicities a 10"^'* — 10~^^ range for ^Be is being 



probed [17[, this would constitute gross overproduction 
of ^Be, which is clearly excluded. 

Figure [2] gives the result of our CBBN calculation at 
Tg < 0.12. Besides rates (HU), ([12]), and ^ calculated 
in this paper, we include the ^Li CBBN ra te [j| with the 
use of the S*- factor properly calculated in [10| . 
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FIG. 2: Logio of the elemental abundances in CBBN for two 
choices of {Yx, tx} input. The top figure represents very long life- 
time Tx and abundance of Yx = 5 X 10~^. This option is excluded 
because ^Be is produced with O(10~^^ — 10"^") abundance. The 
increase in (*BeX~) at low T is due to continuing molecular forma- 
tion. The lower plot represents tx = 2000s and initial abundance 
of Yx = 0.1, suggested by the solution to the ^Li problem Q- For 
these parameters ^Li = 1.3 X 10~^^ , and ^Be = 7 X 10"^*. 



The main result of this calculation is summarized as: 



^Be 



10"^^ X 



[Yx{t = 2 X 10*scc)/10"'^] 



(15) 
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Given that in some cases ^Be is detected below xlO 
17| . and typical abundances Yx > O(10~^) are expected 
from annihilation of at the freezeout, this re- 

stricts the lifetime oi to a few thousand seconds, re- 
inforcing the lithium bound 4] . It is also important that 
^Be is far less fragile than ^Li and therefore is unlikely 
to experience a significant reduction of its abundance in 
subsequent evolution after the Big Bang, making predic- 
tion (|15p especially valuable. Close inspection of Fig. [2] 
reveals a very similar behavior for ^Li and ^Be. Dividing 
the two abundances, we eliminate the dependence on Yx, 
obtaining the following relation, 



BCpi-in^ordial/ Lipi-jmordial — (2 5) X 10 



(16) 



^Li results as a "primordial plateau" , and take seriously 
[very tenuous] hints of elevated levels of ^Be at lowest 
metallicities. It is of course possible that all abnormal- 
ities observed in lithium and beryllium abundances will 
find astrophysical explanations having nothing to do with 
the Big Bang [l^, and only more theoretical and obser- 
vational work in this direction will clarify this issue. 

To conclude, metastable negatively charged particles 
are predicted by many particle physics models, includ- 
ing some variants of supersymmetry. These particles, 
should they live in excess of 1000s, trigger CBBN and 
in this paper we show that sizable amounts of ^Be can 
be generated, and calculated with reasonable accuracy 
from first principles. Somewhat reduced ^Li abundance 
0, strongly enhanced primordial values of ^Li and ^Be, 
and ^Be/^Li'-^^fewx 10"'^ constitute a typical "footprint" 
of CBBN. In this light, further observational studies of 
^Be and ^Li at low metallicities find an unexpected and 
very strong motivation from modern particle physics. 



which is valid as long as 2000s. Intriguingly, this 

is exactly what is observed 17, [l^, if we interpret the 
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